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Abstract 

We consider a system composed of a single artificial atom coupled to a cavity mode. The artificial 
atom is biased such that the most dominant relaxation process in the system takes the atom from 
its ground state to its excited state, thus ensuring population inversion. A recent experimental 
manifestation of this situation was achieved using a voltage-biased superconducting charge qubit. 
Even under the condition of 'inverted relaxation', lasing action can be suppressed if the 'relax- 
ation' rate is larger than a certain threshold value. Using simple transition-rate arguments and a 
semiclassical calculation, we derive analytic expressions for the lasing suppression condition and 
the state of the cavity in both the lasing and suppressed-lasing regimes. The results of numerical 
calculations agree very well with the analytically derived results. We start by analyzing a sim- 
plified two-level-atom model, and we then analyze a three-level-atom model that should describe 
accurately the recently realized superconducting artificial- atom laser. 



I. INTRODUCTION 



Superconducting circuits have gained increased interest in recent years, particularly for 
their possible use in quantum information processing 

to the artificial-atom concept, the idea of placing such an atom in contact with a harmonic- 
oscillator circuit element, which serves as a cavity, has attracted a great deal of attention 
Such circuit-QED systems hold promise for studying various quantum-optics phenomena in 
a highly controllable and easily tunable setting, as well as exploring parameter regimes that 
are inaccessible using natural atoms. 

One of the most intriguing and counterintuitive phenomena in the fields of atomic physics 
and quantum optics is lasing 3|. Given the above-mentioned advantages of superconducting 
circuits for studying atomic-physics and quantum-optics phenomena, it is natural to investi- 
gate superconducting implementations of lasing. Indeed, there have been a number of recent 



theoretical proposals 



9. 



10| and 



and experimental demonstrations of lasing 
population inversion H[ in superconducting systems. 

In Ref. {4] a cyclically manipulated artificial atom is constantly driven into its excited 
state, from which it can relax by emitting a photon into the cavity, thus establishing a lasing 
state. In Ref. j^] an atom that is illuminated by an oscillating field with a properly chosen 
frequency emits photons into a low-frequency cavity. Here we analyze a situation that is 
different from both Refs. {4,6], but is closer to the usual picture of lasing with natural atoms. 
Furthermore, the models that we study are closely related to the experiment of Ref. |9[]. It 
should be mentioned here that similar models have been studied in the past in the study of 



single-atom lasing (see e.g. Refs. 12|, [l3j). A similar model was also analyzed in Ref. 



a, 



but that paper explored different parameter regimes and analyzed different aspects of the 
problem from the present paper. 

Using a transition-rate-based calculation, a semiclassical calculation and numerical sim- 
ulations, we analyze the different possible states of the cavity as the system parameters are 
varied. Each one of the analytic calculations has its advantages. The transition-rate-based 
calculation derives in a transparent manner the lasing suppression condition and the state 
of the cavity deep in the lasing and the suppressed-lasing regimes. The semiclassical calcu- 
lation provides a good approximation for the state of the cavity throughout the lasing state, 
but is not suited to analyze the suppressed-lasing regime, where it turns out that the state 



of the cavity takes the form of a thermal state. For clarity, we start by analyzing a simplified 
two-level-atom model, and we later take the same approach to analyze a three-level-atom 
model that describes more accurately the experiment of Ref. 9J. In particular, we comment 
on a possible experimental implementation of the crossover between the lasing and thermal 
regimes with a superconducting artificial-atom laser. 



II. TWO-LEVEL ATOM 



In this section we analyze the simplified model where the atom contains two energy levels 
only. This model provides a good qualitative understanding of the mechanisms at play 
and the resulting phenomena in the experimental setup of interest to us. The qualitative 
understanding developed in this section will also be useful for identifying the importance of 
the different processes in the more realistic model analyzed in Sec. Ill below. 



A. Model 




K 



FIG. 1: (Color online) Schematic diagram of a two-level atom interacting with a cavity mode. The 
coupling strength for the exchange of excitations between the atom and the cavity is g. The atom 
is biased such that it experiences 'inverted relaxation' from the ground to the excited state, with 
rate T. The loss rate of photons out of the cavity is k. 

We consider the simple system composed of a two-level system interacting with a har- 
monic oscillator (which typically is one mode of an electromagnetic cavity). The system 
is shown schematically in Fig. 1. The Hamiltonian of the combined atom-cavity system is 
given by 

H = ^^a z + hu a^a + ga x (a + , (1) 

where u a is the atom's characteristic frequency, uoq is the cavity's natural frequency, g is 
the atom-cavity coupling strength, a x and a z are the usual Pauli matrices operating on the 



atomic state, and a and a) are, respectively, the annihilation and creation operators acting 
on the state of the cavity. We shall describe quantum states using the notation \n a ,n c ), 
where n a = for the atomic ground state and n a = 1 for the excited state, and n c represents 
the number of photons in the cavity. 

In order to have efficient emission of photons from the atom into the cavity, the atom 
and cavity frequencies must be almost equal. For the remainder of this paper, we shall take 
oj a = uj . We also take this frequency to be the largest frequency (or energy) scale in the 
problem. 

The setup is designed such that the atom's bias conditions cause it to 'relax' from the 

n 

ground state to the excited state, with rate T [15j. In this section, the inverted relaxation 
is assumed as part of the theoretical model under consideration; a similar process will be 
derived from first principles in a realistic three-level-atom model in Sec. III. It is this counter- 
intuitive, inverted relaxation that will provide the mechanism for population inversion, which 
plays a crucial role in the realization of the lasing state. As such, one can say that the (usual) 
threshold condition for lasing action is automatically satisfied in this model. Note that we 
are ignoring any weak relaxation process pushing the atom from the excited to the ground 
state, since such a process would not affect the main points we wish to study. Furthermore, 
since the atom's relaxation rate will be taken to be very large, we shall ignore any additional 
atomic dephasing mechanisms. The cavity is taken to possess a decay rate k. 

An alternative description of the above situation concerning the bias conditions would 
be to say that the cavity is in contact with a heat bath that has a very small and positive 
temperature, while the atom is in contact with a heat bath that has a very small and negative 
temperature. It is worth mentionin g he re that a similar approach (with negative effective 
bath temperature) was used in Ref. 14| to describe an amplification process. 



B. Photon emission and loss rates 



In order to determine the state of the cavity for a given set of parameters, we first 
note that the above model contains a mechanism for photon emission into the cavity and a 
mechanism for photon loss from the cavity. We consider these two mechanisms separately. 

The loss rate of photons from the cavity (i.e. the transition rate from the state \n a ,n) to 
the state \n a , n — 1), where n a represents the state of the atom) is given simply by the decay 



rate k multiplied by the number of photons in the cavity n: 



ri oss = UK. (2) 

Obtaining the photon emission rate requires a somewhat more careful analysis. We first 
consider the situation where there are no or few photons in the cavity. The atom's bias 
conditions constantly push it to its excited state. We can therefore assume the atom to be 
initially in the excited state. If the atom is in its excited state and the cavity has n — 1 
photons, the atom-cavity coupling (with matrix element g^/n) induces dynamics between 
the states |1, n — 1) and |0, n). Since T ^> g, the dynamics will take the form of an incoherent 
process described by the transition rate W\i )n -i)^\o ;n ) = 4ng 2 /T. Any population that starts 
to accumulate in the state |0,n) will quickly relax to the state |l,n), because the atom 
is constantly pushed in this direction by its surrounding environment. These two steps 
complete the transition from the state |l,n — 1) to the state |l,n), or in other words, the 
process of adding one photon to the cavity. Since the upward-relaxation process occurs at a 
very large rate, it can be treated as being instantaneous. We therefore find that the photon 
emission rate (i.e. the transition rate from the state \n a ,n — 1) to the state \n a ,n)) is given 
by the rate of the \l,n — 1) — > \0,n) transition, i.e. 

4ng 2 



emission -p 



(3) 



The photon emission rate therefore increases linearly with n for small values of n [19| . Clearly 
this situation cannot persist for large n, since this mechanism is ultimately limited by the 
atom's relaxation rate T. Indeed, when g^/n becomes comparable to or larger than T, the 
\l,n — 1) «-> |0, n) transitions must be treated as coherent oscillations. One can now argue 
that in the limit of very large n, where the system spends half of the time in each one of 
the two states — 1) and |0, n}), the atom has a chance to incoherently relax from its 
ground state into its excited state only half of the time. In this case the photon emission 
rate asymptotically reaches the value T/2, which it cannot exceed. 

The main advantage of the above derivation of the photon emission rate is its simplicity, 
as well as the simplicity of the resulting expressions. A more detailed analysis of the photon 



emission rate for any value of n is possible, assuming that cavity is in a co 
state. This calculation will be carried out in Sec. II. D below (see also 12 



rerent (i.e. lasing) 



16|). 



C. Lasing condition and possible steady states 



Combing the photon emission and loss rates as functions of photon number n, one can 
obtain the probability distribution of photon number states in the cavity In particular, if 
this probability distribution has a peak for some value of n, the peak value can be obtained by 
locating the intersection point between the emission and loss rates. Some relevant examples 
of such a peak-finding calculation are depicted schematically in Fig. 2. 




FIG. 2: (Color online) Schematic plot of the photon emission rate r em i ss i on (green solid line) and 
the photon loss rate T\ oss (blue dashed line) as functions of the photon number in the cavity n. 
The intersection point (red circle) determines the value of n in the photon-number probability 
distribution with the highest probability. Going from (a) to (d), T is increased, while k is kept 
fixed. 

We treat the atom's relaxation rate T as the tunable parameter, keeping g and k fixed. 
From Eq. ([3]) one can see that small values of T correspond to a large initial slope of the 
emission rate (at n — 0), even though the emission rate reaches the saturation level for a 
relatively small value of n. Figure 2(a) represents this situation. If we increase T, the initial 
slope of the emission rate decreases, but it eventually reaches a larger value. By comparing 
Figs. 2(a) and 2(b), one can see that the peak value of the photon number in the cavity 
increases with increasing pumping rate. This result agrees with intuitive expectations. 

A change of behaviour occurs when T reaches a certain regime, an expression for which 



will be given shortly. As can be seen from Fig. 2(c), the peak value of n starts decreasing 
with increasing V and vanishes at a certain value of V . Beyond this point the value of n 
with maximum occupation probability remains at zero. 

This suppression of lasing action by strong pumping is quite counterintuitive. It can 
be understood in terms of the quantum Zeno effect; at a certain point, the decoherence 
associated with pumping becomes the most dominant effect and inhibits the emission of 
photons from the atom into the cavity [l?], Q|. In the following we shall analyze this effect, 
as well as the state of the cavity in the different regimes, more quantitatively. 

By combining Eqs. (j2J) and (jBD, we find that if 

> *, (4) 

the photon emission rate is larger than the photon loss rate, assuming a small photon number 
in the cavity. Starting with a small photon number, the number increases exponentially in 
time [The growth in photon number continues until the peak value of n is reached, as 
represented by the circles in Figs. 2(a,b)]. If, on the other hand, Eq. (j4]) is not satisfied, the 
loss rate will be higher than the emission rate, and lasing would not occur. Equation 
can therefore be considered a second threshold condition for lasing in this setup. Note that 
population inversion is guaranteed in this model and that all emission from the atom goes 
into the cavity. 

We now consider the situation where the lasing condition (Eq. H]) is satisfied, and we 
analyze the probability distribution of the photon number in the cavity. Deep in the lasing 
regime, we can assume that the emission rate is well approximated by T/2. The loss rate 
is still given by Eq. (j2J). The peak in the photon-number probability distribution therefore 
occurs at 

"max 71 ' (^) 

Note that this steady-state photon number is independent of the atom-cavity coupling 
strength. It is also worth mentioning that this relation remains valid even if T is smaller 
than g. 

The width of the probability distribution can be calculated as follows. The 'probability 
current' from the (n — l)-photon state to the n-photon state is given by 

Wn-^n = (6) 



whereas the probability current in the opposite direction is given by 



W n ^ n -i = riKP n 

r 

= -P n + (n- n max ) KP n . (7) 

Here P n is the probability of having n photons in the cavity. Using the detailed balance 
equation, i.e. W n -i^> n = W n _+ n _i, the above two equations can be combined to give 

Pn Pn—1 Tl Tlmax 



P n 



(8) 



which can be integrated to give the probability distribution 



Pn = P m axeX P {- (n rW)2 |. (9) 

The width of the probability distribution is therefore of the order of ^/n max , as would be 
expected for the lasing state. 

We now turn to the situation where lasing is suppressed, i.e. when 4g 2 < Yk. In the 
linear regime (i.e. when n is small), we can write simple detailed balance equations for the 
probabilities P n ; 

Pn+l _ remission (jl) _ ^9 f~\r\\ 

~Pn~ ~ r loss (n + i) ~ T^' 1 j 

This equation can be identified as the detailed-balance equation for a cavity in thermal 
equilibrium at effective temperature 

(11) 



rri *J 

-LeS 



Tk 



l0 Si —2 



k B 

Note that here we are neglecting the small ambient temperature of the cavity. Using the 
Bose-distribution formula, we find that the average number of photons at the above effective 
temperature is given by 

<»> = ($-i)"- (12) 

Therefore, if we start from large values of T and we gradually decrease it, the average 
number of photons in the cavity starts increasing. This number follows a 1/ rr-type function 
that diverges at the threshold condition. The nonlinearity in the emission rate [see Fig. 2(c)] 
prevents the photon number from diverging at the critical value of T; instead the system 
changes behaviour and enters the lasing regime. 



D. Semiclassical derivation 



In this subsection, we briefly review a mean-field approximation that can be used to find 
an analytic expression for the number of photons in the cavity in the lasing state. We follow 
closely the calculation of Ref. [12| (see also Ref. [181]): we write equations of motion for the 
expectation values of the relevant operators (in the rotating frame for simplicity), and from 
the stationary steady-state solution we extract the number of photons in the cavity. 

We start with the Lindblad master equation for the model under consideration (see 
e.g. Ref. [l4]) 



dp 
~dt 



k apa 



l„ t „ 1 „ t „ 
— a'ap — —pa ] a 



(13) 



where p is the total system's density matrix, and o± are the atom's raising and lowering 
operators. We can now multiply this equation on the left by any operator A and take the 
trace over the density matrix. The result is an equation of motion for the average value of 
the operator A, denoted (A). 

The relevant equations of motion are: 

d (a) 



K 



dt 
d{aJ) 

dt 
d{a z ) 
dt 



r 

9 (a<Tz) ~ g (*-) 

-2g (at <r.. + (UT + ) + r (!-(*,)) 



(14) 



Using the mean-field approximation (i.e. setting (aa z ) = (a) (a z ) etc.), choosing (a) to be 
real and setting the left-hand sides to zero (for the steady-state solution), we find for the 
average number of photons in the cavity (using the relation (n) = (a 2 )) 



(n) 



r 

2k 



El 

v. 



(15) 



This expression is the mean-field approximation of the number of photons in the cavity. It 
predicts that deep in the lasing regime, i.e. when the second term inside the parentheses can 
be neglected, the number of photons will be given by V /2k. It also predicts that the photon 
number will start decreasing with increasing V and will vanish when TK/(Ag 2 ) = 1. Both 
of these results agree with the results of Sec. II. C. Note that the semiclassical calculation 



deals with average values, and therefore is not suited to describe the thermal state (which 
requires knowledge of the probability distribution of the photon occupation number in the 
cavity) . 

E. Numerical calculations 

We now solve Eq. ( Tl3|) numerically for different values of T, keeping g and k fixed. As 
representative quantities that manifest the differences between the lasing and suppressed- 
lasing regimes, we plot in Fig. 3 the average photon number in the cavity (n) and the photon 
number with maximum probability n max as functions of the parameter Tn/(Ag 2 ). 
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FIG. 3: (Color online) Average photon number (n) (blue solid line) and maximum-probability 
photon number n max (green solid line) in the cavity as functions of the parameter Tn/^g 2 ). Note 
that n max corresponds to red circles in Fig. 2. The values g/coo = 8 x 1CP 3 and k/u>o = 5 x 10~ 3 / (2tt) 
were used in the numerical calculations. The red dashed line shows the predictions of Eq. ([5]) in 
the lasing regime, and the red dotted line shows the predictions of Eq. (|12|) in the thermal regime. 
The green line agrees very well with the predictions of the semiclassical calculation (Eq. [15]) . 

The average photon number (n) agrees with the analytic expressions of Sec. II. C [Eqs. §5§ 
and (TI2|) ] away from the threshold on both the lasing and thermal sides. The maximum- 
probability number n max agrees with the quadratic function derived in Sec. II. D (see also 
Ref. throughout the lasing regime. n max coincides with (n) deep in the lasing regime, 
but it decreases faster as the threshold is approached and clearly exhibits an abrupt change 
of behaviour when the threshold condition is crossed. 




FIG. 4: Occupation probability as a function of photon number in the cavity for a point in the 
lasing regime (a; Tn/^g 2 ) = 0.5) and one in the thermal regime (b; TK/(4g 2 ) = 1.25). The system 
parameters are given in Fig. 3. The curve in (a) is fitted very well by a gaussian function, and the 
curve in (b) is fitted very well by a Boltzmann thermal-distribution function. 

In Fig. 4 we plot the probability distribution of the photon number in the cavity for two 
points in Fig. 3, one in the lasing state and one in the thermal state. Apart from a small 
regime around the lasing-suppression threshold, the probability distribution is fitted very 
well by a gaussian function in the lasing regime and by an exponential (i.e. Boltzmann- 
distribution) function in the thermal regime. 



III. THREE-LEVEL ATOM 

We now consider a model that corresponds more closely to the experiment of Ref. Q], 
i.e. a Cooper-pair box coupled to a harmonic-oscillator circuit element. In the analogy 
with conventional lasers using natural atoms, the Cooper-pair box plays the role of the 
atom, whereas the linear circuit element plays the role of the cavity. We follow the methods 
explained in Sec. II above and apply them to this more realistic model with a three-level 
atom. It should be noted here that there has been some work in the past on single-atom 
lasers using three-level atoms or ions The model we consider here, however, provides 

a more accurate description of the experimental situation of main interest to us [9( . 

A. Model 

The Hamiltonian of the system is now given by 

H = — — - (cos + sin 9a x ) + hu a^a + g a z (a + , (16) 



where, as before, oo a is the atom's characteristic frequency and lo is the cavity's natural 
frequency, and these two frequencies are taken to be equal. The angle 9 represents the 
deviation of the atom's bias point from the so-called degeneracy point, and g is the atom- 
cavity coupling strength. In the three-level-atom model, the Pauli matrices a x and d z operate 
on the two active atomic states, with no need to include the third (inert) state explicitly 
in the Hamiltonian (in particular, the energy of the third state does not affect the results 
below). 

A schematic diagram of the photon emission mechanism is shown in Fig. 5, including 
the dissipative processes. The state with N + 1 Cooper pairs in the box can relax to the 
state with iV Cooper pairs and an unpaired electron in the box and an electron added to 
the drain electrode (this transition occurs with rate Ti). This state can relax further when 
the unpaired electron tunnels from the box into the drain electrode (with rate I^). Once 
the box has N Cooper pairs and no unpaired electrons, its coupling to the source electrode 
allows a new Cooper pair to tunnel from the source electrode to the box. 

One can now see how the inverted relaxation process occurs. Using an applied gate 
voltage to the Cooper-pair box, the system is biased such that the state with N +1 Cooper 
pairs in the box is lower in energy than the state with N Cooper pairs, assuming a fixed 
number of electrons in the drain electrode. Under this condition, the state of the box with 
N + 1 Cooper pairs can, on a qualitative level, be identified as the ground state |0), and the 
state with N Cooper pairs in the box can be identified as the excited state |1). When the 
box starts in the state with N + 1 Cooper pairs and one electron tunnels out of the box into 
the drain electrode, the artificial atom goes from the ground state |0) to a third state that 
contains N Cooper pairs and one unpaired electron in the box. Barring coincidences, the 
extra unpaired electron in the box acts as an additional gate voltage, moving the Cooper- 
pair box away from the degeneracy point and from resonance with the cavity. As a result, 
this third level will not be involved in any coherent dynamics and can, for our purposes, be 
considered completely inert. Once in the inert state, the unpaired electron can tunnel from 
the box to the drain electrode and the atom relaxes to its excited state |1), thus completing 
the relaxation process. The fact that the states of the source and drain electrodes change 
during the relaxation processes, which is needed in order to ensure that energy is always 
lowered in each step of the relaxation process, does not need to be explicitly taken into 
account once we have established the mechanism for the inverted relaxation process in the 
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FIG. 5: The different processes involved in lasing for the experiment of Ref. [9|]. The first, second 
and third quantum numbers represent, respectively, the number of Cooper pairs in the box, the 
number of unpaired electrons in the box and the number of electrons in the drain electrode. The 
box is resonantly coupled to a cavity mode, and the two can exchange excitations. The states with 
iV and N + 1 Cooper pairs in the box are coupled because the box is biased in the vicinity of 
the so-called degeneracy point, so that Cooper pairs can tunnel coherently between the box and 
the source electrode. The system's total energy is lowered every time an electron tunnels (in a 
dissipative process) from the box to the drain electrode. The state with N + 1 Cooper pairs and 
a single unpaired electron in the box does not participate in the lasing mechanism. The inset in 
the bottom-left corner of the figure shows the same processes as in the main part of the figure, 
but in the energy eigenbasis of the Cooper-pair box. Because of the mixing between the Cooper- 
pair number states in the energy eigenstates of the box, the different relaxation rates shown in 
that inset obey the relations 70-+2 = Ti cos 2 (#/2), 71^2 = Ii sin 2 (0/2), 72^1 = r2Cos 2 (6'/2) and 
72^0 = r*2 sin 2 (0/2). The inset in the top-right corner of the figure shows a truncated model of a 
three-level atom where relaxation processes take the atom from the ground state (|0)) to an inert 
state (|2)) and then to the excited state (|1)). 



Cooper-pair box. Note that there are additional relaxation processes in Fig. 5, which occur 
because the ground and excited states of the box are superpositions of the states with N 
and N + 1 Cooper pairs in the box. 



B. Photon emission and loss rates 



The derivation of the emission rate is less straightforward in this case than the in Sec. II. B. 
Nevertheless, it can be done, and we carry it out here. 

First we consider the small-n limit, where photon emission can be treated as an incoherent 
process. When the atom is in the state |1), the process by which it emits a photon into 
the cavity still occurs with the rate given by Eq. ([3]), noting that the relevant relaxation 
rate here is Ti (which is the relaxation rate out of the space of active atomic states and 
represents the decoherence rate in that space) and now we have g = go sin 9. One difference 
between the present case and that of Sec. II. B is that now the atom can undergo recurring 
transitions between the states |0), |1) and |2) even without the emission or absorption of 
photons. Therefore we need to include some additional arguments in order to take the 
above fact into account. When the atom is in the state |0), the process by which it absorbs 
a photon from the cavity occurs with the same rate as the one for photon emission. The net 
photon emission rate is therefore given by 

Ang 2 



emission 



[Pi ~ Po) , (17) 



where Pj is the occupation probability of atomic state j. Using the relaxation rates shown 
in Fig. 5, we can derive the probabilities of the different states: 
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P? = e , ~^ 9TT7' (|S) 

which gives 

Pi-Po = • (19) 

cos 2 #+(i + ^)sin 2 # 1 ' 

With this expression we find that the net emission rate in the small-n limit is given by 

r . . = ^ ( ^ \ (20) 

1 emission -ml „ /i \ n I ' V w / 

Deep in the lasing regime, strong coupling with the cavity causes the atom to quickly 
reach equal populations of the states |0) and |1) every time it enters the space of active 



states. Relaxation from the active space to the inert state therefore occurs with rate T\/2. 
Relaxation from the inert state back to the active space occurs with rate IV The atom's 
resetting rate (or Cooper-pair current) is therefore given by (1^/2) x T 2 / (ri/2 + r 2 ), which is 
the rate for a sequence of such recurrent relaxation steps. Since the atom has two possibilities 
when relaxing from the inert state to the active space [state |1) with probability cos 2 (6/2) 
and state |0) with probability sin 2 (#/2)], the net photon emission rate will be given by 

r emlran . <Ii§^ x cos,. ,21) 

1 l/Z + 1 2 

In the following subsection, we use the above rates to describe the state of the cavity. 



C. Lasing condition and possible steady states 



We can now follow the arguments of Sec. II. C with the expressions just derived in 
Sec. III.B and describe different properties of the system. 

Setting the photon emission rate in the small-n limit (Eq. [201) equal to the photon loss 
rate (Eq. [2]), we find the threshold condition 

T 1 k cos 6 

V = cos 2 #+(i + ^) sin 2 9 ( } 

Deep in the lasing regime, equating the emission rate (Eq. [5Tj) to the loss rate (Eq. [2]) gives 
an average photon number of 

i (r x /2) x r 2 



(n) 



cos 6. 



k iV2 + r 2 

In the suppressed-lasing regime, we find a thermal state with effective temperature 

r lK cos 2 #+ (i + ^)sin 2 fl' 



(23) 
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(24) 
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Except for the above modified expressions, the qualitative physical description of the system 
remains essentially the same as the one given in Sec. II. C. 



D. Semiclassical calculation 



We now follow the semiclassical approach to derive the average photon number in the 
lasing state for the three-level-atom model. Using Eqs. (fT4"|) as a template and Fig. 5 as a 
guide for the relevant dissipative processes, we write the equations of motion: 



d (a) 

dt 
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dt 

d(Pi-Pp) 
dt 

d(Pi + Pp) 
dt 
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(25) 



where the 7s are the different relaxation rates, and cr_ transforms the state |1) into the state 
|0). Using the relations 7o_> 2 = ricos 2 (#/2), 71^2 = risin 2 (fl/2), 7 2 _>i = r 2 cos 2 (#/2) and 
72^0 = r2 sin 2 (6 l /2), the steady-state solution of the above equations gives 



1 

2k 



70-^2 - 7l-»2 + 2 72- 
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2->l 



'7^2 7l-+2* 



70^2 + 71^2 + 272^0 + 272^1 



4 



4 



2~2 + 72-»0 + 72^1 



r\ -2r 2 f r 2 /tcos^ 



£1 

2k 



ri + 2r 2 
1 



8g 2 



+ r 2 cos^- 



1 ^ 2r 2 



1 - ^ 

cosfl- I 1 + ^cos 2 fl 

1 + ^ 



(70^2 + 7l->2) 2 K 

4 2gi 

8^ 2 



+ r 2 cos # 



+ 72^1 - 72-*0 



(26) 



Equation fl26|) is clearly more complicated than Eq. (ITB1) . reflecting the more complicated 
nature of the two-step relaxation process in the three-level-atom model. 

If we go deep into the lasing regime, i.e. by neglecting the terms containing Tin/g 2 in 
Eq. ([26"]) . we recover Eq. (1231) . If we equate (n) to zero, we recover the threshold condition 
in Eq. ( [22]) . If we take the case where I\ = I^, Eq. ( 1261) reduces to 



(n) 



£1 

2k 



3 8# 2 V 3 



and the threshold condition is given by 



4 cos 9 



4g 2 3 + cos 2 fl' 
If we take the limit r 2 ^> 17, Eq. ( 1261) reduces to 



(n) 



£1 
2k 



r k 

COS - At f 1 + cos 2 
8g 2 V 



and the threshold condition is given by 

Tk 2 cos 
V ~ 1 + cos 2 0' 



(27) 



(28) 



(29) 



(30) 



For the parameters quoted in Ref. [9j], i.e. g = (2vr) x 44 MHz, T 1 = T 2 = (2vr) x 600 MHz, 
K = (2tt) x 1.3 MHz, 9 = 0.187T, one finds that the ratio T^/^g 2 ) « 0.1 (and (n) « 70), 
with the threshold occuring at rx/t/^g 2 ) « 0.9. This set of parameters is therefore well 
inside the lasing regime. By reducing g and increasing k (e.g. during fabrication), however, 
the boundary between the two regimes seems to be easily reachable. Since the pumping 
rate V is somewhat controllable in experiment, it should be possible to study the transition 
between the two regimes on a single sample. 



E. Numerical calculations 



We solve the quantum-optical master equation relevant to this model [which follows 
straightforwardly from Eq. fflBl and Fig. 5] numerically for different values of Ti, keeping 
g, k and T2/T1 fixed. We plot in Fig. 6 the average photon number in the cavity (n) and 
the photon number with maximum probability n max as functions of TiK,/(Ag 2 ). The main 
features in Fig. 6 are similar to those in Fig. 3, which is an indication that a good intuitive 
understanding of the system can be obtained from the simplified two-level-atom model. The 
curves in Fig. 6 also agree with the analytic expressions given in this section. We do not 
plot the probability distributions here because they look very similar to the ones shown in 
Fig. 4. 



IV. CONCLUSION 



We have analyzed the lasing behaviour of a single artificial atom in a cavity, in particular 
in connection with recent experiments on superconducting charge qubits. Although increased 
pumping strength initially results in a larger photon population in the cavity, increasing 
the pumping rate beyond a certain point starts to suppress the number of photons in the 
lasing state. When the pumping rate reaches a critical threshold value, lasing action is 
completely lost and a thermal state of the cavity is formed. We have analyzed the properties 
of both the lasing and suppressed-lasing (thermal) states. We have used a transition-rate- 
based approach, semiclassical calculations and numerical simulations in our analysis, and 
all three methods give consistent results. Our analysis and results are very relevant to 

n 

the experimentally achieved situation of Ref. |9|, suggesting that experimental tests of the 
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FIG. 6: (Color online) Average photon number (n) (solid lines) and maximum-probability photon 
number n max (dashed lines) in the cavity as functions of the parameter Tin/(4:g 2 ). The values 
g/ojQ = 8 x 10 -4 , k/u>q = 5 x 10~ 4 /(27r) and 9 = ir/3 were used in the numerical calculations. 
The ratio T2/T1 is 1 for the blue (black) lines and 10 for the green (gray) lines. All the numerical 
results agree well with theoretical predictions. The small difference between the solid and dashed 
lines deep in the lasing regime is due to the discreteness of n max . 

phenomena studied in this paper should be possible in the near future. 
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